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RELAXATION METHODS EOR HYPERBOLIC PDE 
MIXED-INTEGER OPTIMAL CONTROL PROBLEMS 

FALK M. HANTEt 


Abstract. We extend the convergence analysis for methods solving PDE- 
constrained optimal control problems containing both discrete and continuous 
control decisions based on relaxation and rounding strategies to the class of 
first order semilinear hyperbolic systems in one space dimension. The results 
are obtained by novel a-priori estimates for the size of the relaxation gap based 
on the characteristic fiow, fixed-point arguments and particular regularity the¬ 
ory for such mixed-integer control problems. As an application we consider a 
relaxation model for optimal flux switching control in conservation laws moti¬ 
vated by traffic fiow problems. 


1. Introduction 

Decision taking, e. g., specifying on/off for some fixed actuator on a plant, is 
a very elementary control mechanism. Critical infrastructure systems such as gas 
pipeline networks, water canal networks or highway traffic networks are for instance 
controlled by switching valves [23], weirs [20] or speed limit signs [15], respectively. 
These decision are possibly to be taken in combination with determining additional 
continuously variable parameters such as the outlet pressure of a compressor in the 
gas network example. Further modeling aspects of such problems are discussed in 
[13]. We will refer to the open-loop optimization of such heterogeneous controls as 
a mixed-integer optimal control problem. In the literature, such problems are also 
called diserete-continuous, hybrid or switching optimal control problems. 

A method for solving such problems numerically should yield integer feasible so¬ 
lutions based on sufficiently accurate approximations of the plant’s dynamics. Due 
to the combinatorial complexity this is a difficult problem in general. In particu¬ 
lar, it becomes a real challenge when the plant is modeled using partial differential 
equations (PDEs) as it is expedient for example in the mentioned infrastructure 
systems in nonstationary scenarios on large scale networks. Thus, several different 
solution approaches have already been proposed in the literature. They can be clas¬ 
sified as total-discretization-, reformulation- and direct relaxation-based and each 
of them have there assets and drawbacks. 

Total discretization of the underlying dynamical system obviously leads to mixed- 
integer nonlinear programs (MINLPs). These become typically large, in particular 
in the PDE case. Hence, solving such MINLPs requires structure exploiting algo¬ 
rithms. If they are available, such methods can provide global optimal solutions. 
This approach is followed for example in [9] for gas networks or in [11] for traffic 
flow. 

Reformulations turn the mixed-integer optimal control problem into—on some 
level—equivalent formulations which can then be solved with existing numerical 
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methods. For example, a variable time transformation method yields a continuous 
formulation [21, 10]. This approach is yet limited to problems governed by ordi¬ 
nary differential equations (ODEs). Optimal switching controls can in some cases 
be sophisticatedly characterized using the viscosity solution of a Hamilton-Jacobi- 
Bellman system, see [7] for ODE problems and [26] for a generalization to PDEs. 
However, in particular in the PDE case, the Hamilton-Jacobi-Bellman system is 
again typically large and numerically a difficult problem. Complementarity-based 
reformulations are used in [3] and have been applied to gas network optimization 
in [4]. The resulting system then require special nonlinear solvers, or additional 
relaxation techniques. 

Direct relaxation of the decision variables obviously turns the problem into a fully 
continuous setting. The resulting problem can then be solved with well-established 
nonlinear optimal control techniques and integer solutions can be recovered from 
rounding. For ODE problems, it has been shown that the relaxation gap, i. e., the 
distance of the optimal value for the relaxed problem and the one corresponding to 
the rounded solution, can be made arbitrary small using relaxation for a convexi- 
fied problem together with a constructive rounding strategy [24]. The so obtained 
epsilon-optimal controls may involve frequent switching, but suboptimal solutions 
can still be obtained including combinatorial constraints limiting for example the 
number of switches [25]. The approach has also been generalized to dealing with 
vanishing constraints [18] and to mixed-integer control problems for abstract semi- 
linear evolutions on Banach spaces using semigroup theory [14]. The generalization 
to Banach spaces covers semilinear PDE cases, but the analysis in [14] assumes a 
certain degree of smoothness of the solution. Due to the inherent discontinuities of 
the integer control interfering with the solution, these regularity assumptions are 
only known to be valid in parabolic cases due to the natural smoothing properties 
of ellipic operators (examples are discussed [14]). In particular, these regularity as¬ 
sumptions are known to fail for classical solutions and are in general very difficult 
to be verified for weak solutions in hyperbolic cases. Unfortunately, all examples 
mentioned at the beginning naturally involve hyperbolic PDEs which are at present 
not supported by the available theory for this approach. 

The main contribution of this paper is to extend the analysis for the direct re¬ 
laxation approach to semilinear hyperbolic PDEs with distributed mixed-integer 
control. Instead of employing semigroup theory on classical Sobolev spaces as in 
[14], we obtain novel a-priori estimates on the size of the relaxation gap using the 
method of characteristics and particular regularity results on semi-classical Sobolev 
spaces. This is an important step to establish relaxation techniques for optimiza¬ 
tion of the mentioned infrastructure systems. Eor instance, transient gas network 
operation can often be accurately modeled using semilinear Euler-equations [I] . For 
an application of the results obtained in this paper to gas-network optimization, 
valve switching is then to be modeled as a distributed control on very small pipe 
sections. This is intended as future work. In this paper, we will discuss the applica¬ 
tion of the method to a relaxation model of nonlinear conservation laws motivated 
by traffic flow control. 

The paper is organized as follows. In Section 2 we give a detailed problem for¬ 
mulation for a semilinear hyperbolic mixed-integer optimal control problem and 
relate it to a relaxed and convexified problem. In Section 3 we estimate the gap 
made by this approach in terms of the integrated difference of two controls. The 
latter quantity is not a norm on the respective control spaces which make the 
estimates technically difficult. The estimates yield a convergence result for the 
relaxation method in Proposition 8. In Section 4, we employ the relaxation ap¬ 
proach to the already mentioned example motivated by traffic flow, where we use 
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an adjoint-equation based gradient-decent algorithm to compute solutions to a re¬ 
laxed problem and we apply rounding strategies in order to verify the predicted 
convergence numerically. In Section 5, we draw a brief conclusion. 


2. Problem Formulation and the Relaxation Approach 


For some real constants L, T > 0, a natural number n, a diagonal matrix function 
A = diag(Ai,..., A„) with Xi : [0, L] ^ R, i G {1,..., n}, normed vector spaces U 
and V and a nonlinear function /: R" x U x V ^ R", we consider a controlled 
system of semilinear hyperbolic PDFs in two variables t (time) and x (a single space 
variable) 

yt + lX{x)y:^ = f{y,u{t),v(t)) on (0, T) x (0, L) (1) 

for an unkown vector function y = (?/i,..., y„)^: (0,T) x (0,L) —>■ R" and two 
controls u: [0, T] —>■ C/ and v: [0, T\ V. We assume that, for some r G {1,..., n}, 
Ai < 0 for j G {1,..., r} and Xi > 0 for z = {r-|-l,..., n} and set y~ = {yi, ..., yr)^, 
= (yr+i, ■ ■ ■ ,yn)^, so that y = {y~,y~^)^. Further, we consider (1) subject to 
boundary conditions 


fy it,L)\ 


/G__ G_+\ {y-{t,0)\ 

{g+. G++) \y+it,L)) 



on (0,T) 


( 2 ) 


and an initial condition 

y{0,x)=y{x) on[0,L], (3) 

where we assume that G _G R’’^’’, G_^ G G^ _G G++ G 

d": [0,T] ^ R’', d+: [0,r] R("-''), and y: [0,L] R". The 

diagonal form in (1) is without loss of generality, because it can be achieved by a 
change of variables for many physical systems, see e. g., [8, 6]. 

Our interest will be to find controls u and v (piecewise smooth and piecewise 
constant, respectively) that minimize a cost criterion 

J{y{T,-))= [ g{y{T,x))dx, (4) 

Jo 

where g: R” ^ R is a non-linear function, subject to the ordinary control constraint 
that 

u{t)€U ( 5 ) 

for some closed (and otherwise arbitrary) set U C U and the discrete control con¬ 
straint 

v{t)€{v\v^,...,v^} ( 6 ) 

for some given values G F, M G N. 

To this end, we note that for any piecewise smooth u and any piecewise constant 
V, the equation (1) is equivalent to 

M 

yt+Ay^ = '^ajf{y,u,v^) on (0,T) x (0, L) (7) 

i=i 

with new controls = aj(t) G {0,1}, j = 1,..., M, satisfying = 1 for 

almost every t. We will in the following estimate the gap (measured in terms of the 
norm of y) of weak solutions for (7) (hence solutions to the original problem) and 
the relaxed problem 


M 

yt + ^yx = '^(3jf{y,u,v^) on (0,r) x (0,L) 
i=i 


( 8 ) 


with controls fdj = j3j(t) G [0,1], j = 1,..., M, satisfying = 1 for almost 

every t. Note that we only replaced the condition “new controls in {0,1}” with the 
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condition “new controls in [0,1]”. Further, observe that the minimization of (4) 
subject to the relaxed problem (8) with the boundary data (2) and (3) and control 
constraints 

M 

(5), /3j(t) e [0,1], j = and '^/3j{t) = l, tG(0,T)a.e. (9) 

i=i 

is an optimal control problem of a seinilinear hyperbolic system in standard form 
which can be solved efficiently with existing numerical methods [16, 2]. 

The main motivation for the a-priori estimates derived in Section 3 is the fol¬ 
lowing observation proved in [24]. 

Lemma 1. Let j3'. [0, T] — >■ [0,1]^ be a measurable function such that — 

1 for almost every t. Then, there exists a piecewise constant function a: [0,T] —>■ 
{0,1}'^ satisfying ^ € [0,T] such that 

max sup 

where At is the length of the largest subinterval where a is taken constant. 

Remark 2. The proof in [24] for the existence of a as in the previous Lemma is 
constructive, giving rise to a numerical method that was also explicitly used in 
[14]. Moreover, it has recently been show that the right hand side in (10) can been 
improved to 0(log(M))At [19]. 


/ P3i.s)-aj{s)ds 
Jo 


<{M - l)At. 


( 10 ) 


We are therefore interested in estimates on the relaxation gap in terms of the 
left hand side of (10). Note that this quantity is not a norm and that these esti¬ 
mates therefore are neither classical nor obvious. For our analysis, we will consider 
solutions of (1) or (8), respectively, subject to the boundary conditions (2) and 
the initial condition (3) in the sense of the forward characteristic flow. Letting 
r^T = [0,r] X [0,L] and Si(t;T, ct), i £ {l,...,n} denote the characteristic curves 
defined by 

4-Si = Aj(si), Si(r;T, cr) = (t,ct), (r, ct) G (H) 

at 

the forward characteristic flow of (8) is obtained as a fixed point y. ^It —t R" of 
the following integral transformation 


^(2/)(A cr) = ('0l(y)(A cr), . . . , 'f’u{y){T, cr)) 
with 'ipi(y){T, cr) defined recursively as 


( 12 ) 


IVl pT 

'0i(y)(Acr) = F*('i/';r,cr) / Pj{t)f{y{t,Si{t;T,<T)),u{t),v^)dt, (13) 

4=1 


where t* = t*{T,a) denotes the intersection time of the curve Si{-;T,a) with the 
boundary of LIt backward in time and Yi(gf;T,<7) is defined as the f-th component 
of 


Vg- g_+\ /v'-(2/)(i:,o)\ /d-(t:)\ 

< l^G+_ G++) \f3+iy){t*,L)) \d+it*)) 


if t* > 0, 


(14) 


.y{t*^s,{t*]T,(j)) ift-=o, 

where 'if~ = [ifi,... ,'ifrY ^ncl = {ifr+i, ■ ■ ■ ■ The so defined solution 

coincides with the usual weak solution of (1) or (8), respectively, subject to the 
boundary conditions (2) and the initial condition (3) in appropriate spaces [12]. 
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3. A-priori Estimates on the Relaxation Gap 

The subsequent analysis is based on a particular regularity result for hyperbolic 
initial boundary value problems obtained in [22] , giving a sufficient condition for the 
solution defined by the forward characteristic flow having bounded distributional 
derivatives along almost every characteristic curve if the initial and boundary data 
is piecewise smooth. To be more specific, for any p,q, ^ S N and any family of 
disjoint open sets Q.m C m = 1,..., /r, we will denote by 

(g)C°(f]™;R«) (15) 

m—1 m—1 

the set of functions h defined on the closure of Um=i with image in R'? so that 
their restriction to ilm belongs to the classical Sobolev space R'^) or the 

Banach space R'^), respectively. 

Then, we make the following assumptions. 

Hypothesis 3. The components of K, Xi : [0, L] —>• R, i € {1,..., n} are Lipschitz- 
continuous. The functions y —>• f{y,u,v^) are smooth and satisfy f{0,u,v^) = 0 
for all j = 1,..., M and u € U. The function u —>■ f{y, u, v^) is locally Lipschitz- 
continuous on U for all j = 1,...,M and y € Y. Further, for all T > 0, 
there exist finitely many points 0 = tq < ti < ... < tk-i < tk = T and 
0 = xo < xi < ... < Xv-i < x,^ = L so that the initial and boundary data 
satisfies y G cci+i;R"), d~ G lT^’^(ri,Ti+i;R’'), and d'^ G 

These assumptions yield the following wellposedness result. 

Lemma 4. Under Hypothesis 3 and for any piecewise smooth control it: [0, oo) U 
there exist constants T, K > 0 such that (8) subject to the boundary conditions (2) 
and the initial condition (3) admits a unique solution y in (^([O, T]; L°°(0, L; R”) n 
L^(0,L;R")) for all piecewise smooth controls fi: [0,T] —>■ [0,1]'^. In particular, 
this solution is given as the (unique) fixed point of (12)-(14) resulting from a strict 
contraction with contraction constant ^ on y = C([0,T];L^(0,L;R”)) equipped 
with the norm 

lll/llt = sup / \yi{t,x)\dx. (16) 

i G [0,T] 2=1 ^ 

Proof. Hypothesis 3 implies that the initial and boundary data, respectively, satis- 
fiesy G L^’(0,L;R"), d- G LP{0,T;W) and d+ G LP(0,T; R^""’')) for allp G {l,oo} 
and all T > 0 and that / is locally Lipschitz-continuous in y. Noting that all feasible 
fd are uniformly bounded by one in the sup-norm, the result follows from classical 
fixed-point arguments [8], [6]. □ 

In order to obtain further regularity properties of the solution, we adapt the fol¬ 
lowing definitions from [22]. For any piecewise smooth control fd'. [0,oo) —>■ [0,1]'^ 
with discontinuities at {0i}fdi, let E^{(d) denote the union of all forward charac¬ 
teristic curves Si generated by A (and their reflections at the boundaries) which lie 
in rioo = [0, oo) X [0, L], which emerge from the boundary points 

{(0,xo),...,(0,a:,.), (tq, 0),..., (ri^, 0), (tq, A),..., (rjf, 0)} (17) 

and their intersection points with the sets {9i} x [0, L] for all z = 1,..., oo. Further, 
for fc G N, let E^{fd) be the union of E^~^{fd) and all forward characteristic curves 
(and their boundary reflections) emerging from intersection points of characteristics 
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that defined Let Et{P) be the closure of all points in U^i^^(/3) l~l 

Supposing that 

Et{P) n {t} X [0,L] is nowhere dense in {t} x [0, L], t G (0,T], (18) 

Et{P) is defined by finitely many discrete curves, which divide VLt up into finitely 
many simply connected open sets fl™, m = 1,..., /r, // S N. We set 

\ Erm := (g) C°(L!t \ Erm := (g) C\n^). (19) 

m—1 m—1 

With this notation, we note the following additional regularity of the L^-solution 
for piecewise smooth data and controls. 

Lemma 5. Under Hypothesis 3 and for any piecewise smooth control u: [0, oo) -G U 
and any piecewise smooth control fd: [0, oo) ^ [0,1]"^, there exist a constant T > 0 
(less or equal to the constant obtained in Lemma 4) such that (18) holds for the set 
Et{/3) defined above. Moreover, it holds that the solution y of (8) subject to the 
boundary conditions (2) and the initial condition (3) obtained in Lemma 4 satisfies 

y e W(’\nT \ Erm n CliLlr \ Et^P))- (20) 

Proof. The result follows directly from [22, Theorem 3.1 and Remark 3.3] applied 
iteratively on time intervals where d~, d"*", u and (3 are jointly continuous. □ 

Moreover, we will need the following fixed point argument. 

Lemma 6. Consider a Banach space {X, || • ||x). Let At := {x € X : ||a;||f < oo} 


for some norm || • ||| satisfying 

c||a;||x < Ikllt < c||a;||x, x G X. (21) 

Let (j), Ip: X ^ X be continuous mappings such that 

||(/)(x)-(/)(C)||t < ^lk-Cllt> IIV'(a::)-V'(C)llt < ^l| 2 ;-Cllt> x,CgX. ( 22 ) 

Suppose that there exists a constant C > 0 so that 

mx*)-^P{x*)\\x<C, (23) 

where x* G X denotes the (unique) fixed-point of f in X. For the (unique) fixed 

point G X of tp in X it then holds 

\\x*-C\\x<2<r^cC. (24) 

Proof. By the triangular inequality and (21) we have 

||X* - Cllt < mx*) - f^{x*)\\x + Mx*) - ^(Ollf (25) 

Using (22) and (23), this yields 

lk*-Cllt<cC'+i||x*-C|lt, (26) 

and equivalently (24) using again (21). □ 


With the above auxilary results, we can now prove the following a-priori estimate 
on the difference of two solutions corresponding to two different controls with a 
bounded integerated difference. 

Theorem 7. Assume Hypothesis 3 and let u : [0, oo) U, (3: [0, oo) [0,1]^ and 
p3'. [0, oo) —>■ [0,1]*^ be piecewise smooth controls. Moreover, assume that for some 
given e > 0 and T sufficiently small, 

Pjis) - fij{s)ds 



max sup 


< e. 


(27) 
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Then, there exists a constant C = C{u, P) > 0 (independent of P) so that 

\\y{u,P)-y{u,P)\\Y <C£, (28) 

where y(u,P) and y[u,P) denote the -solutions of the relaxed problem (8) subject 
to the boundary conditions (2), initial condition (3) and the control constraints (9) 
with controls {u, P) and {u,P), respectively. 

Proof. Let T be less or equal to the constant T obtained in Lemma 5 for the 
fixed control {u,P). Moreover, let ip denote the integral transformation (12)-(14) 
associated to the fixed control {u, P) and let p denote the corresponding integral 
transformation associated to the fixed control {u,P). By Lemma 4, there exists 
K > 0 such that both transformations are strict contractions with contraction 
constant ^ on with respect to the norm || • ll-j. defined in (16), possessing the 
unique fixed points y = y{u, P) and y = y(u, P), respectively. 

We will now show existence of a constant C such that 


/ l'0i(y)(T,cr) - P^iy)iT,(T)\da <Ce, i = l,...,n, t G [0,T] 

Jo 

holds for the fixed point y = y(u,P). 

We have for r € [0, T] and 0 < = t* (r, a) < t 

[ \'Pi{y){T,(j) - pi{y){T,a)\da <\\Yi{tjj;T,a)-Yi{p;T,a)\\ + 


M L 


■ 


da. 


da < 


[Pj{iJ) - Pj{'d)]f{y{'d, Si{iJ; r, a)), m(i?), u*) dd 
Integration by parts for the integral in d yields 

[/3j(i?) - Pj{d)]f{y{d, s^{'d;T,a)),u{d),v^)dd 

Dsf{y{'&, Slid-, T, a)), uid), v^) \ ddda 
\fiyiT,SiiT-,T,a)),uiT),v^)\da + 
\fiyit*,Siit* ; T, a)),u{t*),v^)\da. 


M 

U 

M I.L 

E 


1=1 
M „L 

u 

fr 


0 






By assumption (27), this estimate becomes 

[Pjid) - Pjid)]fiyid,Siid-,T, a)), u{d), ) dd 


M .L 

Yf 


1=1 • 


da < 


M 

^E 

1=1 


CT pL 


'tt ^0 


\D^f{y{'d,Siid-,T,a)),u{d),v^)\da + 


IfiyiT, Si(T;r, a)),uiT),v^)\ + |/(?/(t■, Si(l■; r, a)),uit*),v^)\da 


(29) 


(30) 


(31) 


(32) 
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The chain rule yields 


n pL 


E 

k^l' 


Dykf Siiid; r, a)), uid), ) 


Dufiviid, Si(29; T, C 7 )),u(d),v’) 


D^u{'d) 


D^yk{'d,Si{'d-,T, a)) 
da. 


+ (33) 


The assumptions in Hypothesis 3 on /, the assumed regularity of u and the reg¬ 
ularity result on y = y{u,fi) in Theorem 5 now yields that all terms in absolute 
values appearing in the right hand side of (33) are essentially bounded. Since all 
remaining terms in the right hand side of (32) are finite and independent of j3, there 
exists a constant C independent of /3 such that 

E/ / - ^jW]f(y('^^Si{-d\T,a)),u(§),v^)d'd 


da < Ce. 


(34) 


Moreover, for \\Yi{'ip]T,a) — li((/); t,ct)||, we can use (14). For i S {l,...,r} and 
(r, a) such that t* > 0 for all cr G [0, L], we obtain 


\\y^{'^;T,a) - Y,{(j);T,a)\\ < HGH 


(2/)(i3‘>0) - (/> (y)(t*,0)\ 

(y) (tt , 0) - <(>+ (y) (t* , 0 )J 


(35) 


and a similar estimate is obtained in the other cases. We can then repeat the above 
estimates for each of the components of the vectors in the right hand side of (35). 
Since t* < t and since the estimate becomes trivial for t* = 0, the existence of a 
constant C independent of P such that (29) follows from induction. 

Finally, we note that X = G([0,T];L^(0,L;]R")) is a Banach space for which 
(21) holds with c = and c = 1. So the desired estimate (28) follows (29) and 

Lemma 6 with X = y and C = 2c“^cnC' = 2e^'^nC. □ 


A combination of Theorem 7 and Lemma 1 yields a method for solving the mixed- 
integer optimal control problem stated in Section 2 subject to Hypothesis 3 up to 
any requested accuracy. We state this method in form of the following proposition 
and remarks. 


Proposition 8. Assume Hypothesis 3 and suppose that for some sufficiently small 
T the minimization of (4) subject to the relaxed problem (8) with the boundary 
data (2) and (3) and control constraints (9) admits a piecewise smooth optimal 
control {u, /3) with associate optimal value J*. Let be a sequence of controls 
: [0,T] {0,1} given by Lemma 1 obtained from /? for a sequence At^ —>■ 0 as 

k ^ oo. Define : [0,r] V by v^{t) = oi^{t)v^. Then y(u,v^) converges 
strongly to y{u, (3) as k ^ oo. In particular, for norm-continuous cost functions J, 

lim J{y{u,v^)) = J*. (36) 

k—fco 

For Lipschitz-continuous J, the convergence speed is at least linear in At^. 

Proof. This follows from applying Lemma 1 with a = and Theorem 7 with 
P = and e = (M — l)At*^ for each fixed k. □ 

Remark 9. We may relax the assumption of existence of a piecewise smooth opti¬ 
mal control (u, /3) in Proposition 8 by existence of a piecewise smooth suboptimal 
control (m, 0) with J* = J{y{u, /?)). The statement of the Proposition then does not 
change otherwise, but it clearly yields a slightly weaker conclusion. We may also in¬ 
clude combinatorial constraints on the integer controls, e. g., bounding the number 
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of switches between the modes. Then, we do in general not obtain the convergence 
(36), but the relaxation method then still yields asymptotically a monotonic de¬ 
creasing sequence J{y{u,v^)) of suboptimal solutions as /c oo. Moreover, the 
relaxation method can also deal with state constraints. Details of these extensions 
are discussed in [14]. 

Remark 10. Note that Theorem 7, Proposition 8 and Remark 9 are stated for 
sufficiently small T > 0 guaranteeing the existence of the solution for the control 
(m,/ 3) both in the sense of Lemma 4 and Lemma 5. However, this limitation can 
be removed under additional assumptions on /, A and the system’s dimension. 
For example, is well-known that if the mappings y i—>■ f{y, m, v^) and x i—>■ A(x) are 
Lipschitz-continuous, then the L^-solution (coinciding with the L°°-solution) in the 
sense of Lemma 4 exists for all T > 0. Further, in the important case of n = 2, we 
have Et{/3) = and hence Et{I3) satisfies (18) for all T > 0. In that cases, 

for example. Theorem 7 and Proposition 8 hold for arbitrary T. Nevertheless, for 
n > 3, even if the L^-solution exists and coincides with the L°°-solution globally in 
time, Et{P) may violate (18). For n = 3, an example is also given in [22]. Hence, 
in general, for n > 3, (18) must be explicitly checked. 

4. Application to flux switching control of conservation laws 
We consider the integer controlled semilinear hyperbolic system 

Vt+^x = 0 , 

(37) 

Ct + a'^Vx =we {1,2} 

with 5(17,1) = 5(5,2) = At > 0 and 0 ? such that — rf > 0. In 

characteristic variables 51 =5-1-0^ and 52 = the system (37) can be written 

as (1) with A = diag(a, —a) and a nonlinear function /(y, u, v) independent of u. 
For sufficiently small k, the system (37) is an approximation of the control system 

Vt ±Wx = 0 (38) 

where the control just consists of switching the sign in the flux function of the 
conservation law (38). The approximation holds in the sense that, for fixed v, up 
to second order in k 

Vt±^vl = K{{a^-r]'^)r]^)^, (39) 

see, [17, 5]. Such flux switching control problems appear for example when traf¬ 
fic flow modeled by conservation laws is supposed to be optimized by switching 
dynamic speed limit signs as in [15]. Since Burger’s equation is a typical test prob¬ 
lem for traffic flow, this example shall demonstrate that the relaxation method 
investigated in this paper is a very efficient way to solve such problems with the 
advantage of being supported by a convergence theory and being applicable up to 
the same discretization levels that can be handled by optimal control techniques 
for hyperbolic systems without integer confinements. 

For our example, we consider the initial data 5(0, x) = r]o{x), ^(0, x) = ^5(50,1)+ 
^5(50,2) = 0, for a given 50, periodic boundary conditions = r]{t,L), 

^(t, 0) = ^(t, L) at the end points of the interval [0, L] and consider the minimization 
of a tracking type cost functional 

•^(5) = ^^ \r]{T,x)-f]{x)\^dx (40) 

for a given reference solution fj. Assuming 50 and 5 being piecewise on (0,L) 
(which is not restrictive in most applications), we obtain from Lemma 5 that the 
L^-solution of (37) is piecewise on (0, L) for each t G [0,T]. Since in this 
example we have n = 2, (18) is satisfied for all T > 0, cf. Remark 10. Hence, for 
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any T > 0 such that the usual L^-solution exists, r]{T, ■) and fj are, by Sobolev- 
embeddings, also in L^(0, L) and the evaluation of J is therefore well-posed for any 
piecewise constant control v. [0,T]—>-{1,2}. 

According to the relaxation method presented at the beginning of Section 3 
motivating Proposition 8, we introduce two new controls 01,02: [0,T'] —>■ {0,1}, 
such that oi(t) -|- 02(<) = 1, and consider their relaxation /3i,/32: [0,T'] —)• [0,1] 
subject to /3i(t) + /32{t) = 1. Using the latter constraint yields /?2(i) = 1 — /3i(i), 
so that the the relaxed control problem (8) contains only a single control function 
P{t) = /3i(t). Moreover, using that we have/3(t)(^—5(77, l))-|-(l—/3(t))(^—(7(77, 2)) = 
-I- the relaxed problem (8) reads for this example 

= 0 , 

^t + a^'n^ = -\y^-£.)■ 

The optimization problem that determines j3 in Proposition 8 is 
min J(77) subject to 
T =0; 

6+a^?7x = K"^(/3(t)77^ - - C) (42) 

77(0,2:) = 770(x), ^( 0 ,a;)= 0 , r]{t,0) = r]{t,L), L) 

m e [0,1]. 

We solve (42) using an adjoint-equation based gradient decent method. To this 
end, letting (p, q) being the solution of the following adjoint equations 

- pt - a^qx = K ~^ q { 2 P ( t ) - 1)77, 

-qt-Px = -K~^q (43) 

p(T, x) = -ir]iT, x) - ri{x)), q{T, x) = 0, p(t, L) = p{t, 0), q(t, L) = q{t, 0), 
the derivative of the reduced cost function J(/3) = J{r]{l3)) can be obtained as 

J'W) = -f -b(??, 2)-5(77,1)] dx = / -p^dx. (44) 

Jo ^ Jo ^ 

As in [17, 2], we use a first order finite-volume in space and implicit Euler in time 
splitting scheme for the discretization of (41) and (43). The integrals in (40) and 
(44) are evaluated using the trapezoidal rule. 

As a test problem, we consider L = 2tt, T = 3, Pq(x) = 2x(i x € [0, L], 

with x-y denoting the characteristic function of the set 7 in [0,i], and fj{x) = 
1 — sin(x), X G [0,L]. We discretize the space domain [0,L] by Nj; = 300 cells, 
set a = 5, K = l.Oe—08 and choose a CFL consistent time discretization step size 
with the CFL constant 1/2 for the time interval [0,T]. The discretization level 
corresponds to a MINLP with 4.276.800 unknown real and 7.138 unkown binary 
variables. 

With the above adjoint-based approach, we could easily computed a piecewise 
constant /3* with J(/3*) = 0.086 up to a first order optimality of l.Oe—08 using an 
interior point method. The numerical results corresponding to Proposition 8 are 
reported in Table 1. The computed controls, the initial state and the corresponding 
final time states are plotted in Figure 1 and 2. 

In this example, we observe numerically that the relaxed problem exhibits very 
nice bang-bang structure. This structure is captured by the rounding strategy of 
Lemma 1 for At = 0.25 or smaller. We have observed similar results for other 
initial data and control targets. Hence, for the application to flux switching control 
of the form (37), the relaxation method shows even better convergence properties 
than those predicted in Proposition 8. In particular, in this example, we do not 
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k 

At^ 

J{v^) 



1 

1.00 

0.467 

0.381 

4.44 

2 

0.50 

0.396 

0.310 

3.62 

3 

0.25 

0.086 

0.000 

0.00 

4 

0.125 

0.086 

0.000 

0.00 

5 

0.0075 

0.086 

0.000 

0.00 


Table 1. Numerical results for the test problem on flux switching 
control for conservation laws discussed in Section 4. 



Figure 1 . Optimal control for the relaxed problem and the integer 
control obtained by sum up rounding strategies for different At for 
the test problem in Section 4. 


observe frequent switching in the optimal integer control. This motivates even fur¬ 
ther investigation of the structure of solutions to hyperbolic mixed-integer optimal 
control problems. 


5 . Conclusions 

Our analysis and numerical results show that certain PDE mixed-integer optimal 
control problems of hyperbolic type can be solved successfully und very efficiently 
using methods based on relaxation and rounding strategies. 
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